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Abstract 



This paper presents observations of a protoplanetary disk around Herbig Ae star HD163296 in 12 CO ( J=l-0), 
12 CO ( J=3-2), 13 CO ( J=l-0), and 13 CO (J=3-2) emission lines. Double-peaked emission profiles originating 
PLh ■ from the rotating circumstellar disk were detected in all the lines. The disk parameters were estimated from model 

' calculation in which the radial distribution of temperature or surface density inside the disk has a power-law form. 

The surface density should be sufficiently high so that the disk is optically thick for all the CO lines, as discussed 
in previous studies based on interferometric observations . The temperature and outer radius of the disk were also 
confirmed to be consistent with the previous results. Taking advantage of difference in position of the photosphere 
among the CO lines, we revealed temperature distribution in vertical direction. The temperature of 12 CO ( J=3-2) 
(— > ' emitting region is about twice higher than that of any other CO emitting region; the former is 58.5 ± 9.5 K while the 

latter is 31 ± 15 K at 100 AU from the central star, suggesting that there are at least two distinct temperature regions. 
The best fit temperature for 13 CO (J=l-0) that should trace the deepest region of the disk is even lower, implying 
that there is also a different temperature region in deep inside of the disk. Such vertical temperature distribution in a 
disk was identified both in T Tauri and Herbig Ae stars (e.g., DM Tau, AB Aur, and HD31648), and this should be a 
common feature in protoplanetary disks. 
■ Keywords: stars: individual(HD 163296) — stars:planetary systems: protoplanetary disks — stars: pre-main- 

' sequence — radio lines: stars 

JO ' 1. Introduction Temperature distribution along the vertical direction of the 

\ disk should also be taken into account when one examines how 

Herbig Ae (HAe) stars are widely recognized as interme- the gas and dust inside the disk evolve. Theoretical studies 

diate mass counterparts of low mass T Tauri stars (TTSs) predict that flared circumstellar disks should have two or more 

(Waters & Waelkens 1998). In the last few decades, it has layers with distinct temperature in vertical direction that are 

been revealed that these stars are commonly accompanied by formed by diffusion of stellar radiation scattered at the disk sur- 

*_h . a circumstellar disk like TTSs by direct imaging at infrared face (Chiang & Goldreich 1997; D'Alessio et al. 1999; Inoue, 

and radio wavelengths (Fukagawa et al. 2004; Wisniewski et Oka, & Nakamoto 2009). To examine vertical structure of the 

al. 2008; Mannings & Sargent 1997; Mannings, Koerner, & disks based on observations, it requires examples as many as 

Sargent 1997; Pictu, Dutrey, & Kahane 2003). These disks are possible, though several studies have already tried to make 

believed to be the precursors of planetary systems around HAe vertical temperature or density distributions unveiled (Panic 

stars, and are the important targets to study the formation pro- et al. 2008; Dartois, Dutrey, & Guilloteau 2003). HD163296 

cess of a planetary system around an intermediate-mass star is one of the best targets for this purpose because it exhibits 

and to construct a general scenario for the formation of plane- strong emission from the disk in many CO isotopologue lines, 

tary systems. Recent advanced instruments have revealed their HD163296 is a star with A3V in spectral type (Meeus et al. 

physical properties, such as the disk mass, the radial temper- 2001). Its luminosity, mass, age, and distance are 30 Lq, 

ature and surface density distributions, the outer radius of the 2.3 A/©, 4 Myr, and 122 pc, respectively (van den Ancker, 

disk, and the velocity field. The number of well-studied disks de Winter, & Tjin A Djie 1998). Observations at infrared 

around HAe stars, however, is still insufficient to establish a and radio wavelengths revealed the existence of a circumstellar 

firm understanding of their physical conditions and their statis- disk with its mass of 0.028 M@ and 45° in inclination angle 

tics; only a few of limited examples of their structure have (Mannings & Sargent 1997; Isella et al. 2007; Hughes et al. 

deeply been analyzed (e.g., Panic et al. 2008; Pictu, Dutrey, 2008). HD163296 is categorized in group II by the feature of 

& Guilloteau 2007). spectral energy distribution (SED) and expected to have a self- 
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shadowed disk (Acke & van den Ancker 2004; Meeus et al. 
2001; Dullemond & Dominik 2004). More recently, imaging 
studies with IRAM/PBI, SMA, and VLA both in continuum 
and in the 12 CO, 13 CO, and C ls O emission lines were made 
(Isella et al. 2007), arguing that the disk radius is 550 AU and 
that the gas kinematics can be well explained by Keplerian ro- 
tation. 

This paper presents multi-line observations of HD 163296 
with the Nobeyama 45-meter (NRO 45 m) radio telescope 
and the Atacama Submillimeter Telescope Experiment ( ASTE) 
sub-millimeter telescope, and discusses vertical temperature 
distribution inside of the disk. The outline is as follows: Details 
of the observations are described in section 2, and the results 
are shown in section 3. In section 4, physical properties of 
the disk are estimated by model fitting in which power-law 
distributions of temperature and surface density are assumed. 
Interpretation of the model fitting results as well as the compar- 
isons with previous interferometric observations are discussed 
in section 5. 

2. Observations 

2.1. 12 CO ( J =1-0) Mapping Observations 

Observations of 12 CO (J=l-0) emission line were per- 
formed in winter 2006 by the NRO 45 m telescope, operated 
by Nobeyama Radio Observatory (NRO), a branch of National 
Astronomical Observatory of Japan (NAOJ). To separate the 
emission originated from the disk associated with the cen- 
tral star from surrounding components such as an envelope 
or a remnant cloud, a profile map including the stellar posi- 
tion and its vicinity with sufficient spacial resolution is deci- 
sive in identifying the origin of the emission. To obtain such a 
map, the multi-beam receiver, 25 Beam Array Receiver System 
(BEARS: Sunada et al. 2000; Yamaguchi et al. 2000) that is 
capable of simultaneously detecting 25 different locations in 
double sideband (DSB) operation was employed as the fron- 
tend and tuned at the rest frequency of 12 CO ( J=l-0), 1 15.271 
GHz. The half-power beam width (HPBW) is 15" at this fre- 
quency, corresponding to 1800 AU at a distance of 122 pc 
from the Sun. The beam grid spacing of BEARS is 40" (« 
5000 AU) in both R.A. and Dec. directions, but this is much 
larger than the HPBW and therefore is insufficient to isolate the 
star+disk system from the ambient. To fill these gaps among 
the beams, we also took data at positions shifted by 20" in both 
R.A. and Dec. directions from the central star. The integration 
time for the observation was 7 hours. We only show 80" x 
80" (« 0.05 pc x 0.05 pc) region near the star that covers a 
typical size of the envelope. A digital spectrometer with 32 
MHz in bandwidth and 37.8 kHz in frequency resolution was 
used as the backend. Correction for atmospheric absorption 
during observations was made by the chopper-wheel method 
(Kutner & Ulich 1981), and the intensities were obtained in 
antenna temperature, T£, in kelvin. Correction for the main 
beam efficiency, ?7 m b, was made when we compared model 
calculation results (section 4) with the brightness temperature, 
Tmb = Tl/n mh , where r) mb = 0.39 for 12 CO (J=l-0) emis- 
sion line. The DSB system temperature in 12 CO ( J=l-0) map- 
ping observation was between 250 K and 600 K. The accuracy 
for telescope pointing was regularly checked against the SiO 



maser VX Sgr., which was about 30' away from HD163296, 
and pointing deviation was achieved in less than ±3" during 
whole observation period. Data reduction and analysis were 
made with the NEWSTAR software package developed by 
NRO, which is a front end of AIPS developed in the National 
Radio Astronomy Observatory (NRAO). 

Observations were made in position-switching mode. It was 
found that the extended CO (J=l-0) emission was present 
around the region of interest, hence careful search is needed 
to find an off point as close to HD163296 as possible. The 
search for off points was conducted with BEARS based on the 
fact that the systemic velocity of HD163296 is 6 km s" 1 (Dent, 
Greaves, & Coulson 2005; Isella et al. 2007). The emission 
free position applicable as an off point was finally found at 
about 10' north from HD 163296. Figure 1 is a profile map to- 
ward this position, centered at 17 h 56 m 20?0 and -21°46 m 23?2 
(J2000) and covering 80" x 80" (« 0.05 pc x 0.05 pc) region, 
showing that there is no prominent emission detected around 6 
km s" 1 in U LS r. 

2.2. 13 CO ( J=l-0) Emission Line Observations 

Observations of 13 CO (J=l-0) emission line were con- 
ducted by the NRO 45 m telescope in winter 2008. As will 
be shown in section 3.1, the 12 CO (J=l-0) emission de- 
tected around 6 km s _1 in Vlsr at the central star is isolated. 
Therefore, the 13 CO (J=l-0) observations were made only at 
the stellar position by a pair of single-beam receivers that en- 
able to simultaneously detect both the polarization components 
in single sideband (SSB) mode. The receiver was tuned at 
13 CO ( J=l-0) rest frequency of 1 10.201 GHz and the HPBW 
is 15" at this frequency. The ?y m b in 2008 at this frequency was 
0.4. Acousto Optical Spectrometer (AOS) whose bandwidth 
and frequency resolution are 40 MHz and 20 kHz respectively 
was used. The system noise temperature during observation 
was between 300 K and 600 K in SSB mode. The pointing 
was regularly checked by the same manner as in 12 CO ( J=l- 
0) observations described in section 2. 1 . The same off point as 
in the 12 CO ( J=l-0) observations was used. All observational 
data when the wind speed exceeded 5 m s" 1 were flagged out 
before obtaining the final profile of 13 CO (J=l-0) emission. 
The resultant total integration time was 8.5 hours on source. 
Data reduction and analysis were made in the same manner as 
in 12 CO ( J=l-0) observation in section 2.1. 

2.3. 12 CO and 13 CO (J=3-2) Observations with ASTE 

Observations of 12 CO and 13 CO (J=3-2) line were car- 
ried out in July 2006 with Atacama Submillimeter Telescope 
Experiment (ASTE), a 10 m sub-millimeter telescope in Chile 
(Ezawa et al. 2004). HPBWs of the ASTE telescope at 345.79 
and 330.59 GHz were 22" and 23", respectively, which cor- 
responds to about 2700 - 2800 AU at the distance to the tar- 
get, and the ?y m b was 0.6 around these frequencies. SIS mixer 
receiver having DSB response was employed. Typical atmo- 
spheric opacity at 220 GHz toward the zenith was 0.05 dur- 
ing the observations, and the system noise temperature of the 
telescope was 250 — 300 K in DSB during observations. The 
digital auto-correlator with 1024 frequency channels was con- 
figured so that the total bandwidth was 128 MHz, resulting 
in a frequency resolution of 125 kHz, or ~ 0.11 km s _1 at 
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Fig. 1. Profile map of the off positions by the 12 CO ( J=l-0) observation. The central position is (a, 5)=(17 h 56 m 20?0, 
the angular distance between adjacent points is 40" in both R.A. and Dec. directions. 



-21°46 m 23?2) in J2000. and 



the observed frequencies. Telescope pointing was checked 
every 1.5 — 2 hours by cross scan of Jupiter, and the error 
was proved to be < 3". All the spectra were obtained in a 
position-switching mode: the location of the off-position was 
(Aq, AS) = (-20', +20') from the star. In 12 CO {J =3-2) 
observations, we took spectra not only at the stellar position 
but also at the four adjacent points, 22" apart from the stellar 
position in R.A. or Dec. direction, to check out whether the 
emission of the stellar position is really isolated. 13 CO (J=3- 
2) data were taken only toward the stellar position since the 
12 CO (J=3-2) at the systemic velocity (Vlsr ~ 6 km s _1 ) was 
proven to be isolated, as will be shown in section 3.2. Data 
reduction and analysis were conducted in the same manner as 
in 12 CO (J=l-0) observations described in section 2.1. 



3. Results 

3.1. 12 CO(J=l-0)and 13 CO(J=l-0) 

Figure 2 shows a 12 CO ( J=l-0) profile map within the 80" 
x 80" region centered at HD 163296. A clear double-peaked 
profile with 4.8er is seen only at the center panel at 6 km s _1 in 
Vlsr that matches with the systemic velocity of HD 163296 
(Dent, Greaves, & Coulson 2005; Isella et al. 2007). The 
double-peaked profile has peak intensity of 1 .0 K in T m \, and 
3.2±0. 1 km s _1 in FWHM velocity width. The achieved mean 
rms noise level after data reduction was 0.21 K. Figure 3a dis- 
plays the same profile on the central panel in figure 2 in an 
enlarged view. 

Several positions in figure 2, especially 40" away to the west 
from the center, show an emission feature at Vlsr = 6 km s _1 , 
but there seems no emission above the 3er level at this veloc- 
ity at four positions adjacent to the star (20" apart from the 
star). To confirm whether there is no significant emission at 
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Fig. 4. (a) Averaged 12 CO (J=KJ) profile of the four locations 20" 
away from the central star in figure 2. (b) The 12 CO (J=l-0) profile 
toward HD163296 in figure 2. (c) A profile after subtracting (a) from 
(b). 

6 km s _1 in these regions, these four profiles were averaged 
into one spectrum (figure 4a). There is no emission above 
the 3cr level at 6 km s _1 in figure 4a. The profile detected 
at the stellar position (figure 4b, or the central panel in figure 
2) is compared with that after subtracting the profile in figure 
4a (figure 4c). Although the noise level in figure 4c increases 
due to subtraction, there is no remarkable difference between 
these two profiles. These comparisons strongly suggest that the 
double-peaked emission component at the star is isolated and 
originates from the compact circumstellar disk in Keplerian 
rotation. The components at Vlsr = 6 km s _1 in south and 
west regions > 40" apart from the star in figure 2 might be the 
emissions from a remnant cloud partially remained around the 
star+disk system. 

The dip at 8 km s _1 around the star+disk system (figure 2) 
is due to weak emission seen at the off points, especially in 
the northern region (figure 1); in position switching mode, the 
emission at the off point appears like an absorption dip. The 
dip dose not affect further analysis since it does not corrupt the 
original emission profile at 6 km s _1 . In addition, the strong 
emission components between 10 and 15 km s _1 and at 20 km 
s _1 appear regardless of any positions: these two components 
probably come from foreground or background clouds that ex- 
ist coincidentally along a line of sight. 

Figure 3b shows the profile of 13 CO ( J=l-0) emission line. 



Again, a double-peaked profile was clearly detected. The peak 
intensity and velocity width are 0.36 K in T m b and 3.8±0.2km 
s _1 at FWHM, respectively. The achieved rms noise level is 
0.061 K. Since the 12 CO ( J=l-0) and 13 CO ( ,7=1-0) emissions 
are detected in the same velocity range, both emissions should 
have the same origin. The intensity ratio between 12 CO( J=l- 
0) and 13 CO( J=l-0) is less than 5, suggesting that the disk is 
optically thick. Isella et al. (2007) evaluated the optical depth 
of the CO emission lines and concluded that the disk is opti- 
cally thick. Our model fitting shown in section 4.2 also pro- 
vided that the disk is optically thick for both the CO lines. 

3.2. 12 CO (J '=3-2) and 13 CO (J '=3-2) 

Figure 5 shows a 12 CO ( J=3-2) profile map within the 44" 
x 44" region centered at HD163296; the central panel of figure 
5 is also shown in figure 3c in an enlarged view. The achieved 
rms level was 0.122 K. A clear double-peaked profile at ~6 
km s _1 is seen only in the center panel. The peak intensity 
and velocity width are 1.35 K in T m b and 3.1±0.1 km s -1 in 
FWHM, respectively. The emission components between 10 
and 15 km s _1 and at ^20 km s" 1 in Vlsr appear regardless 
of positions as in the case of 12 CO ( J=l-0) (figure 2). 

Figure 3d presents the 13 CO ( J=3-2) spectrum at the stellar 
position, showing its double-peaked profile. The peak inten- 
sity and velocity width are 0.38 K in T m b and 3.85 km s _1 in 
FWHM, respectively. As in the case of the ,7=1-0 lines, the 
disk seems optically thick for ,7=3-2 line from the intensity ra- 
tio between 12 CO ( J=3-2) and 13 CO ( J=3-2)(see also section 
4.2). 

4. Model Fitting 

Model fitting is one of the effective methods to estimate 
physical properties of disks around HAe/Be stars (Dullemond, 
Dominik, & Natta 2001; Dullemond 2002; Testi et al. 2003; 
Natta et al. 2004). We fitted all the observational profiles, 
which are two CO isotopologues' emission lines with two dif- 
ferent rotational transitions, with those calculated by the disk 
model (details in section 4.1). We adopted a classical disk 
model which includes power-law forms of temperature and 
surface density distributions in radial direction with inner and 
outer cutoffs (Kitamura et al. 1993). The results are compared 
with the ones obtained by millimeter and sub-millimeter inter- 
ferometers (Isella et al. 2007, details in section 5.1). Although 
the interferometric observations can provide detailed geomet- 
rical structure with high spacial resolution, it might be unable 
to detect the emission from the diffuse gas uniformly extending 
away due to spatial filtering effect. The single dish telescopes 
with large diameter, on the other hand, have higher sensitivity 
to the weak extended emission which interferometers can not 
detect well. Interferometric and single dish observations are, 
therefore, complementary to each other, and comparing their 
results is a fruitful way to gain a better understanding of the 
disk nature. 

4.1. Disk Model 

Several assumptions for the disk structure and parameters 
are applied in the model fitting. First, the disk is assumed 
to have Keplerian rotation as was proven to be the cases in 
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Fig. 2. Profile map of the 12 CO (J=l-0) line around HD163296. The map center corresponds to the stellar position (a, <5)=(17 h 56 m 21?29, 
— 21°57 m 21. s 9) in J2000, and the angular distance between adjacent points is 20" in both R.A. and Dec. directions. 



other HAe stars such as HD169142 and HD31648 (Panic et 
al. 2008; Mannings, Koerner, & Sargent 1997). Isella et al. 
(2007) pointed out the existence of Keplerian rotating disk 
around HD 163296 from "a butterfly shape" seen in its posi- 
tion velocity diagram. Secondly, the temperature distribution 
over the disk, T(r), is assumed to follow a power-law form of 
the radius with its index of —0.5, as described by the following 
equation: 



T(r)=T 



-0.5 



(1) 



100AU 

where r is the radius in cylindrical coordinate system centered 
on the star and To is the temperature at r = 100 AU (Hayashi, 
Nakazawa, & Nakagawa 1985; Beckwith et al. 1990). This as- 
sumption implies that the temperature is constant along the ver- 
tical direction (z) at a certain r. In reality, a circumstellar disk 
should have temperature distribution along the z-direction. An 
elaborate disk model (e.g., Chiang & Goldreich 1997; Tanaka, 
Himeno, & Ida 2005) predicts that the surface regions where 
the stellar radiation can penetrate and heat directly will have 
higher temperature than that in the interior regions, forming 2- 



layer temperature structure. Single layer temperature model, 
however, can still be applicable because an optically thick line 
should have "photosphere" at each r and its characteristic tem- 
perature can uniquely be determined. That is, vertical temper- 
ature distribution can be traced by a series of a different tem- 
perature of each CO emission line. 

Thirdly, the radial distribution of surface density over the 
disk, S(r), is also assumed to follow a power-law form of the 
radius with its index of —1.0 (Isella et al. 2007), as described 
by the following equation: 



£(r) = S n 



(2) 



. 100AU 

where Eo is the surface density at r = 100 AU. Fourthly, we 
assumed that density distribution in vertical direction can be 
achieved by hydrostatic equilibrium. The density distribution 
p(r, z) is, therefore, expressed by 



p(r,z) = p(r,0)exp 



H(r) 



(3) 



where p(r, 0) is the density at the mid-plane of the disk, and 
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Fig. 3. Profiles detected in 12 CO and 13 CO lines toward the stellar position of HD163296. (a) in 12 CO (J=l-0), (b) in 13 CO (J=l-0), (c) in 12 CO 
( J=3-2) and (d) in 13 CO (.7=3-2). 



H(r) is the scale height given by 



H(r) = 



l2r 3 k B T{r) 
Gil/* TO 



(4) 



where &b is Boltzmann constant, G is gravitational constant, 
M* is the mass of the central star, and m is the mean mass of 
gas molecules. Thus, p(r, 0) in equation (3) becomes 

E(r) 



p(r,0) 



(5) 



Fifthly, local thermodynamic equilibrium (LTE) is assumed. 
The number density of H2 at the position of photosphere of 
each CO line can be estimated by the equation (3). As shown 
in section 4.2, the number density of H2 at the 12 CO ( J=3-2) 
photosphere is calculated to be ~ 10 7 cm -3 at r = 100 AU, and 
it is higher than the typical value of critical density for 12 CO 
(J=3-2) (~5x 10 4 cm -3 at 60 K, fitting temperature of this 
line discussed in section 4.2) by more than two orders of mag- 
nitude. LTE assumption is sometimes inadequate for analysis 
of protoplanetary disks since it tends to overestimate the ex- 
citation temperature for a higher transition with a high critical 
density. However, Pavlyuchenkov et al. (2007) evaluated many 
types of model calculations and found that LTE approximation 
seems valid for the cases of J=4-3 or a lower transition of 
CO. The LTE assumption is, therefore, applicable in our model 
calculation for the J=l-0 and J=3-2 transitions of 12 CO and 



13 CO. 

Finally, the abundance ratios of H 2 to 12 CO and 12 CO to 
13 CO are assumed to be 10 4 and 60, respectively, which are 
the typical values in the interstellar medium (Frerking, Langer, 
& Wilson 1982; Wannier, Penzias, & Jenkins 1982). The disk 
is assumed to be heated only by a stellar radiation and accretion 
heating is neglected. This is a reasonable approximation for a 
disk with modest accretion rate such as HD 163296 case (7.6 x 
10~ 8 M Q yr _1 ; Garcia Lopez et al. 2006). 

All parameters regarding the object, listed in table 1, should 
be taken into account in the model calculation for the best accu- 
racy, but some of them rarely affect the disk property. Emission 
profiles are mostly characterized by three parameters, To, Eo, 
and the outer radius of the disk (r out ), and these are regarded 
as variables. The other parameters are fixed to be constants as 
listed in table 1 . 

In actual calculations, Eo was fixed in each run, and the pos- 
sible fitting range of To was searched. Eo was set to be be- 
tween 0.003 and 3 g cm~ 2 . The search for To in upper range 
of Eo is stopped at 3 g cm -2 because To in each CO line is 
well confined due to asymptotic behavior of fitting tempera- 
ture in large Eo region and the result remains almost the same. 
Intensity is affected by both To and Eo in the case of optically 
thin disk and by only To in the case of optically thick disk. 
Among three variables, r out can be guessed from the separa- 
tion between two peaks in the profiles obtained from obser- 
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vation, since the Keplerian velocity at r out (V(r ou t)) roughly 
corresponds to half of the separation between the two peaks of 
the profile. r ou t =1000 AU is adopted as an initial guess and 
then adjusted. 

4.2. Fitting Results 

The best fit temperatures and outer radius at 0. 1, 0.3, 1, and 3 
g cm~ 2 in So where fitting solution for all the CO lines simul- 
taneously exist are provided in table 2, and the fitting results are 
shown in figure 6. Note that the observed profiles shown in fig- 
ure 6 were re-sampled with velocity resolution of 0.08 km s _1 
for 12 CO (J=l-0) line and 0.11 km s _1 for the others to di- 
rectly compare with the model calculations. As shown in table 
2, the fitting temperature of 12 CO ( J=3-2) shows the highest 
value among all the CO lines. It is essential to ascertain the 
goodness of the solution obtained in the model fitting. Figure 7 
shows fitting temperature at several specific So between 0.003 
and 3 g cm -2 . We regarded the results as an acceptable fit 
when the following two conditions were met simultaneously: 
(i) to exclude the case of large deviation in amplitude of pro- 
files, the residual in absolute value integrated over the range 
of T m b > 2<r should be less than the uncertainty (Ict) for the 
integrated intensity, (ii) to exclude the case of large deviation 
in line shape, more than 80% of the channels over T m b > 3cr 
should agree with the model calculations within 2a. 

The best fit T for the 12 CO (J=3-2) line in table 2 are be- 



tween 58.5 ± 9.5 K, and this range is significantly higher than 
those estimated for the other CO lines, 31 ± 15 K when So is 
greater than 0.1 g cm -2 . This difference in best fit To implies 
that the disk is composed of at least two layers with distinct 
temperatures. The uncertainty of fitting temperature is smaller 
in higher surface density or large optical depth in the disk. In 
lower surface density, the fitting temperature range becomes 
drastically larger or no solution exists; 13 CO ( J=l-0) emis- 
sion line, for instance, dose not have any fitting temperature 
at So less than 0.1 g cm~ 2 . The fitting errors in temperature 
at 0.1 g cm~ 2 in S are more than 30 % for 13 CO (J=l-0) 
line and more than 20 % for 13 CO ( J=3-2) line, respectively. 
The ones at greater than 0.3 g cm~ 2 in So, on the other hand, 
are less than 20 % for all the CO lines and the fitting tempera- 
tures are gradually confined. The So that is greater than 0.1 g 
cm -2 obtained by model fitting corresponds to > 10 (~ 5 even 
at r ou t = 550 AU) in optical depth, and therefore the disk is 
surely optically thick for all the CO lines. 

5. Discussion 

5.1. Comparison with observation by interferometers 

High-resolution images of HD 163296 were taken with in- 
terferometers at multi-wavelengths (the IRAM/PBI, SMA, and 
VLA) both in continuum and in the emission lines of 12 CO, 
13 CO, and C ls O (Isella et al. 2007). They compared the ob- 
servational results of the CO emissions with those predicted by 
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Table 1. Parameters in model calculation 



Invariables Variables 





2.3 M Q * 


T 


d 


122 pc* 


So 





4 5 °ti 


r out 


r- m 


0.1 AU 




P 


l.ott 




q 


0.5 




X{ 12 CO) 
X( 12 CO)/X( 13 CO) 


1/10000 




60 





* van den Ancker, de Winter, & Djie (lyyS). 
t Isella et al. (2007) 

* Hughes et al. (2008) 




FLSRtkm S" 1 ] FLSRtkm s-'] FLSR[km s-'] FLSR[km s->] 



Fig. 6. (Upper panels) Comparisons between the line profiles obtained by our observations (solid line) and best-fit model calculations when So = 0.1 g 
cm -2 (green lines). Best-fit values for To and r ou t in the model calculations are listed in table 2. (Lower panels) Difference between the observed and 
calculated profiles shown in the upper panels. 



Table 2. Best Fit Parameters at r out = 550AU 



So [g cm" 2 ] 


T [K] 

12 CO(J=3-2) 13 CO(J=3-2) 12 CO(J=1-0) 13 CO(J=1-0) 


0.1 

0.3 
1 

3 


63+° 5 38lg 35t 3 261" 
59t| 33±| 32^ 21±| 
55t 5 3 30tl 30t 3 3 19±| 
53^ 28^ 29^ 18±| 
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Fig. 7. Best-fit values of Tq with uncertainties for the acceptable fit in model calculations as a function of Sq . The panel (a) shows the cases of 



: 0.003 — 3 g cm , and the panel (b) magnifies the lower right part of (a). 



self-consistent disk model to estimate physical disk properties. 
Although the assumptions in model calculations were slightly 
different from each other, the temperature and minimum sur- 
face density of the disk derived in section 4.2 showed a reason- 
able agreement with the values obtained by Isella et al. (2007); 
the temperature at 100 AU and the minimum surface density 
at 10 AU from the central star were estimated to be 26 lg K 
and 1 g cm -2 from our 13 CO ( J=l-0) observations and they 
were within the possible range of the fitting values obtained by 
Isella et al. (2007): 30±10 K and A-tf g cm~ 2 , respectively. 
The radius derived from single dish observations, which gener- 
ally have higher sensitivity to weaker emission, was 550±100 
AU and was consistent with the one derived by interferometric 
observations. 

5.2. Thermal Structure of the Disk 

As described in section 4.2, the disk is arguably optically 
thick in CO emissions. Taking advantage of different mass 
opacity of isotopologues whose emission lines reflect the tem- 
perature at each photosphere having different depth, the verti- 
cal thermal structure of the disk can be analyzed. As shown 
in figure 7, model fittings in the case of 12 CO ( J=l-0), 13 CO 
( J=l-0), and 13 CO ( J =3-2) emission lines were succeeded in 
almost the same temperature when Eo > 0.1 g cm -2 , but the 
12 CO ( J=3-2) results require almost double of the fitting tem- 
perature of the other three lines. These different temperatures 
can be explained by temperature distribution along the vertical 
direction, as predicted by theoretical studies (e.g., Chiang & 
Goldreich 1997). 

To examine the vertical temperature distribution in detail, we 
derived the column densities of H2 that satisfy optical depth for 
each CO emission line to be unity (N Tv —1) from 



N T =1 



3fc B T cx AU f 



gas 



1 



8n 3 Bfi 2 cos9 (J+l) 



exp 



hBJ(J+l) 



1 — exp 



hv 



(6) 



where B is the rotational constant of a molecule, (1 is the 
permanent dipole moment, T ex is the excitation temperature, 
AV^ as is the velocity width of a gas molecule, is the incli- 
nation angle, h is Planck constant, J is a rotational quantum 
number, and i/q is the frequency of transition (Scoville et al. 
1986). Note that N Tu =1 is the value obtained by integrating 
the number density from the disk surface to the disk interior 
along a line of sight. 

In this paper, we adopted AVj, a s = AVth for simplicity, 
where AVth is the thermal velocity width of gas molecules ex- 
pressed by 



AV gas = AU th = 



l8\n2k B T(r) 
mco 



(7) 



Turbulent velocity (14urb) can also contribute to AVgas. 
Hughes et al. (201 1) derived a turbulent line width of ~ 0.3 km 
s _1 for the disk around HD163296 by fitting profile with high 
spectral resolution. Previous observations also revealed that 
I4 ur b in disks around TTS or HAe star is comparable to AVt n ; 
0.07-0.15 km s" 1 in DM Tau (Dartois, Dutrey, & Guilloteau 
2003; Simon, Dutrey, & Guilloteau 2001) and 0.38 km s" 1 in 
AB Aur (Pictu, Guilloteau, & Dutrey 2005). In such cases, 
however, the line profile calculated by the model is not al- 
tered significantly by AVWb because the line shape is mainly 
determined by Keplerian rotation velocity. This assumption 
probably produces a factor of uncertainty for the estimation of 
N Tu —1 but it does not significantly affect the following discus- 
sion. 

N Tvq= i for each CO emission line provides the information 
about the vertical location of the photosphere of each CO emis- 
sion line in the disk. Figure 8 shows JV T =i and vertical loca- 
tions of photospheres of each CO emission line. The results are 
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provided in the case of To = 30 K and 60 K, which are the rep- 
resentative values of the low temperature interior and the high 
temperature surface, respectively, as described in section 4.2. 
To derive the location of each photosphere, on the other hand, 
we select S = 0. 1 g cm~ 2 as the representative case since this 
value is consistent with the results by Isella et al. (2007). In 
both cases shown by figure 8, N Tu = i results in the descend- 
ing order of 12 CO (J=3-2), 12 CO (J=l-0), 13 CO (J=3-2) and 
13 CO ( J=l-0), implying that the photosphere of 12 CO ( J=3- 
2) emission is located uppermost among the four CO lines. The 
total geometric cross section of grains per hydrogen atom (cm) 
in interstellar medium is ~10 -21 cm 2 at visible wavelength 
(Stahler & Palla 2004), hence the column density that makes 
the optical depth for the stellar radiation unity is 10 21 cm~ 2 . 
The grazing angle is assumed to be 0.05 in Inoue, Oka, & 
Nakamoto (2009) calculation, hence the column density of the 
upper layer should be about 5 x 10 19 cm -2 . This column den- 
sity is consistent with the one above the photosphere of 12 CO 
( J=3-2) at r > 100 AU, as shown in figure 8. N Tvo=1 of 12 CO 
( J=l-0) decreases faster than that of 12 CO ( J=3-2) at r > 300 
AU in the case of To = 30 K because CO molecules is not suffi- 
ciently excited to the J=3 level due to lower temperatures. The 
same phenomenon is seen in 13 CO as well. From the compar- 
isons between the fitting temperature and the vertical locations 
of the photosphere of each CO line, the general trend is that 
the farther away the layer is from the mid-plane of the disk, 
the warmer it reaches. This is quite reasonable because the gas 
around the disk surface should be heated more by the radiation 
from the central star. 

When we derive N Tv —i in figure 8, we assumed that the 
disk temperature is uniform in vertical direction. These con- 
siderations, however, might be invalid as a consequence of ver- 
tical temperature distribution; the 12 CO ( J=3-2) emission line 
originates from an upper layer of the disk, but this could also 
contribute significantly to the other CO lines. We therefore 
evaluate how strongly such an upper layer affects the other CO 
emissions. Figure 9 shows how strongly the 12 CO (J=l-0) 
is emitted from above the photosphere of 12 CO (J=3-2). In 
this figure, we first calculate the optical depth of 12 CO ( J=l- 
0) emission line in the region above the photosphere of 12 CO 
( J=3-2), and then the intensity of 12 CO ( J=l-0) at the fitting 
temperature was compared with that originated from the re- 
gion. The contribution of high temperature layer is less than 
10% of the 12 CO (J=l-0) intensity and can be safely ne- 
glected. Therefore, the estimated vertical location of each CO 
emission layer in figure 8 is validated. 

As shown in figure 7, the temperature of 12 CO (J=l-0) 
and 13 CO (J=3-2) layers are almost the same when So > 0.1 
g cm~ 2 . The temperature of 13 CO (J=l-0) layer seems to 
be even lower. Such vertical thermal structure of irradiated 
accretion disks around TTS has been theoretically predicted. 
D'Alessio et al. (1999) constructed the detailed vertical struc- 
ture models of irradiated accretion disks and extensively ex- 
plored the dependence of the structure and emission properties 
on mass accretion rate, viscosity parameter, and disk radius. 
Inoue, Oka, & Nakamoto (2009) also propose a 3-layer disk 
model. The middle layer can form where the optical depth for 
the radiation reprocessed by the upper layer is less than unity. 



Our results suggest the existence of such a three-layer struc- 
ture. 

6. Summary 

This paper presents observational results of the disk around 
HD163296 in 12 CO ( J=l-0), 12 CO ( J=3-2), 13 CO (J=l-0), 
and 13 CO (J=3-2) emission lines. Double-peaked profiles 
originated from a rotating circumstellar disk were successfully 
detected in all the CO lines. Physical parameters of the disk 
such as temperature distribution, surface density distribution 
and outer radius are estimated by model fitting utilizing a disk 
model. These physical values obtained by single dish observa- 
tion were confirmed to be consistent with the results of inter - 
ferometric observations (Isella et al. 2007). 

The disk must be optically thick for all the CO lines. Taking 
advantage of difference in position of the photosphere among 
the CO lines, we revealed temperature distribution in verti- 
cal direction. It is proven that there are at least two distinct 
temperature layers, possibly three layers, from the model fit- 
ting. The temperature in the uppermost 12 CO ( J=3-2) emit- 
ting layer is estimated to be 58.5 ± 9.5 K at 100 AU from the 
central star, and it is about twice higher than that in the inner re- 
gions emitting the other CO lines, 31 ± 15 K at 100 AU. Since 
the vertical temperature distribution inside of a disk is com- 
monly suggested for both TTS and HAe stars (DM Tau, AB 
Aur, and HD31648), such a temperature distribution may be 
ubiquitous in protoplanetary disks (Pictu, Dutrey, & Guilloteau 
2007). The temperature and density structure of protoplane- 
tary disks are still open questions. Better knowledge such as 
the radial and vertical density distribution will be provided by 
high-resolution observations in upcoming ALMA era. 
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Fig. 8. (upper panels) Column densities in ATfH^) that are required to become t vq = 1 for each CO emission line. The radial temperature dependence 
is set to be T(r) = 30(r/100AU) -0 ' 5 [K] (left) and T(r) = 60(r/100AU) -0 ' 5 [K] (right), respectively. Note that the vertical axises are reversed, 
(lower panels) Vertical distance from the mid-plane of the disk normalized by the scale height of the local temperature in the case of So = 0.1 g cm -2 . 
These vertical locations of each CO line correspond to the location corresponding to N Tij0=1 shown in the upper panel. 




Fig. 9. Influence of the higher temperature layer on the CO ( J=l-0) emission by higher temperature layer. The solid line represents the intensity of 12 CO 
(J=l-0) emission at the fitting temperature and surface density (To = 30 K and So = 0. 1 g cm -2 ). The dashed line represents the intensity of 12 CO 
( J=l-0) emission line from the region where the optical depth of 12 CO ( J=3-2) is less than unity when To = 60 K. 
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